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Abstract 

A proposed high energy Future Circular Hadron-Electron Collider would provide sufficient energy in a clean environment to probe 
di-Higgs production. Using this channel we show that the azimuthal angle correlation between the missing transverse energy and 
the forward jet is a very good probe for the non-standard hhh and hhWW couplings. We give the exclusion limits on these couplings 
as a function of integrated luminosity at a 95% C.L. using the fiducial cross sections. With appropriate error fitting methodology 
we find that the Higgs boson self coupling could be measured to be ^ -O^-o i7(o i2) expected Standard Model value at 

yfs = 3.5(5.0) TeV for an ultimate 10 ab“* of integrated luminosity. 
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1. Introduction 

The 125 GeV particle discovered by the ATLAS and CMS 
experiments [1-5] has been established as a spin-0 Higgs bo¬ 
son rather than a spin-2 particle [6, 7]. The measurements of its 
couplings to fermions and gauge bosons are being updated con¬ 
stantly and the results confirm consistency with the expected 
Standard Model (SM) values [4, 8-10]. However, to establish 
that a scalar doublet does indeed break the electroweak (EW) 
symmetry spontaneously when it acquires a nonzero vacuum 
expectation value, v, requires a direct measurement of the Higgs 
boson self coupling, A. The minimal SM, merely on the basis of 
the economy of fields and interactions, assumes the existence of 
only one physical scalar, h, with 7^"' = 0^^. Although Ref. [8] 
has ruled out a pure pseudoscalar hypothesis at a 99.98% con¬ 
fidence limit (C.L.), the new particle can still have a small CP- 
odd admixture to its couplings. 

Theoretically, the Higgs boson self coupling appears when, 
as a result of electroweak symmetry breaking in the SM, the 
scalar potential y(<l>) gives rise to the Higgs boson self interac¬ 
tions as follows: 

y(0) = -H T((1)' <1))2 ^ ^mlh^ + Avh^ + (1) 

where A = dsM = /(2v^) 0.13 and d) is an 51/(2)/, scalar dou¬ 
blet. Eor a direct and independent measurement of A we need 
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to access double Higgs boson production experimentally. How¬ 
ever, this path is extremely challenging and requires a very high 
integrated luminosity to collect a substantial di-Higgs event 
rate, and an excellent detector with powerful background re¬ 
jection capabilities. On the theoretical side we need to also take 
into account all vertices involved in the process that are sensi¬ 
tive to the presence of new physics beyond the SM (BSM). 

There are various proposals to build new, powerful high en¬ 
ergy e^ef e p and pp colliders in the future. We have based our 
study on a Future Circular Hadron-Electron Collider (ECC-he) 
which employs the 50 TeV proton beam of a proposed 100 km 
circular pp collider (ECC-pp), and electrons from an Energy 
Recovery Linac (ERL) being developed for the Large Hadron 
Electron Collider (LHeC) [11, 12]. The design of the ERL is 
such that the e p and pp colliders operate simultaneously, thus 
optimising the cost and the physics synergies between e p and 
pp processes. Such facilities would be potent Higgs research 
centres, see e.g. Ref. [13]. The LHeC and the ECC-he config¬ 
uration are advantageous with respect to the Large Hadron col¬ 
lider (LHC) (or ECC-pp in general) in terms of (1) initial states 
are asymmetric and hence backward and forward scattering can 
be disentangled, (2) it provides a clean environment with sup¬ 
pressed backgrounds from strong interaction processes and free 
from issues like pile-ups, multiple interactions etc., (3) such 
machines are known for high precision measurements of the 
dynamical properties of the proton allowing simultaneous test 
of EW and QCD effects. A detailed report on the physics and 
detector design concepts can be found in the Ref. [11]. 

The choice of an ERL energy of Eg - 60 to 120 GeV, with an 
available proton energy Ep - 50 (7) TeV, would provide a cen¬ 
tre of mass (c.m.s.) energy of yfs ^ 3.5(1.3) to 5.0(1.8) TeV 
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at the FCC-he (LHeC) using the FCC-pp (LFiC) protons. The 
FCC-he would have sufficient c.m.s. energy to probe the Higgs 
boson self coupling via double Higgs boson production. The 
inclusive Higgs production cross section at the FCC-he is ex¬ 
pected to be about five times larger than at the proposed 100 km 
circular collider (FCC-ee). 

This article is organised as follows: We discuss the process 
to produce the di-Higgs events in an e p collider and the most 
general Lagrangian with all relevant new physics couplings in 
Section 2. In Section 3 all the simulation tools and the kine¬ 
matic cuts that are required to study the sensitivity of the in¬ 
volved couplings are given. Here we also discuss the details of 
the analyses that has gone into the study. In Section 4 there is 
a discussion on the validity of the effective theory considered 
here. And finally we conclude and draw inferences from the 
analysis in Section 5. 

2. Formalism 

In an e p collider environment a double Higgs event can be 
produced through: (1) the charged current process, pe“ ^ 
hhjve, and (2) the neutral current process, pe“ —» hhje~, if 
there are no new physics processes involved. The SM back¬ 
ground will cloud each of the processes greatly, and it will be a 
formidable task to separate signal from backgrounds. Here we 
study the charged current process because the signal strength of 
this is superior to the neutral current process. Hence we show in 
Fig. 1 the Higgs boson pair production, at leading order, due to 
the resonant and non-resonant contributions in charged current 
deep inelastic scattering (CC DIS) at an e p collider. As seen in 
Fig. 1, the di-Higgs production involves hhh, hWW and hhWW 



Figure 1: Leading order diagrams contributing to the process pe —> 
hh jve with q = u,c,d,s and q' = d, s, u, c respectively. 

couplings. Note that the hWW coupling will be extensively 
probed at the LHC, where its value conforms to the value pre¬ 
dicted by the SM [4, 8, 10]. Through vector boson fusion Higgs 
production mode at the LHC, a BSM analyses to determine the 
CP and spin quantum numbers of the Higgs boson has been 
studied in Refs. [14-17]. The authors of Ref. [18] have shown 
the sensitivity of new physics contributions in hWW couplings 
at e~p colliders through a study of the azimuthal angle correla¬ 
tion for single Higgs boson production in pe —^ hjve with an 
excellent signal-to-background ratio based on the h ^ bb de¬ 
cay channel. Since we do not have any direct measurement of 
the Higgs boson self coupling {hhh) and quartic (hhWW) cou¬ 
pling, there can be several possible sources of new physics in 


the scalar sector. This article studies for the proposed FCC- 
he sensitivity of the Higgs boson self coupling around its SM 
value including BSM contributions by considering all possible 
Lorentz structures. In order to make it a complete study we 
also retain the possibilities for hWW couplings that appear in 
the di-Higgs production modes. 

Following Refs. [18, 19] the most general Lagrangian which 
can account for all bosonic couplings relevant for the phenom¬ 
enology of the Higgs boson sector at the FCC-he are the three- 
point and four-point interactions involving at least one Higgs 
boson field. It can be written as: 
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Here g® yi - 1,2, and g(...) are real coefficients corresponding 
to the CP-even and CP-odd couplings respectively (of the hhh, 
hWW and hhWW anomalous vertices), W^y = d^W^ - dyW^ 
and W^v - \e^ypa-WP°'. In Eq. (2) is parametrised with a 
multiplicative constant with respect to Tsm as in Eq. (1). Thus 
the Higgs self coupling A appears as in the expression 

for y(0). Clearly, in the SM = 1 and all other anoma¬ 
lous couplings vanish in Eqs. (2)-(4). The Lorentz structures 
of Eqs. (2)-(4) can be derived from the S U(2)l <S>U{1)y gauge 
invariant dimension-6 operators given in Ref. [19]. 

The complete Lagrangian we work with is as follows: 
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The most general effective vertices take the form: 
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The momenta and indices considered above are of the same or¬ 
der as they appear in the index of the respective vertex T. For 
example, in the vertex Thw-w* the momenta of h, and 
are p\,p 2 and p 3 respectively. Similarly, p 2 and /13 are the 
indices of and W^. Using the above effective field the¬ 
ory (EFT) approach a study has been performed as an example 
for di-Higgs production in vector boson fusion at the LHC in 
Ref. [19]. 

3. Simulation Tools and Analysis 

We begin our probe of the sensitivity of these couplings 
by building a model file for the Lagrangian in Eq. (5) using 
FeynRules [19], and then simulate the charged current dou¬ 
ble Higgs boson production channel pe^ ^ hh jvg (see Eig. 1), 
with h further decaying into a bb pair,* in the ECC-he set up 
with ^ 3.5 TeV. Our analysis starts with optimising the 
SM di-Higgs signal events with respect to all possible back¬ 
grounds from multi-jet events, ZZ-i-jets, hbb - 1 - jets, hZ + jets 
and ffH-jets in charged and neutral current deep-inelastic scat¬ 
tering (DIS) and in photo-production^, taking into account ap¬ 
propriate / 7 -tagged jets and a high performance multipurpose An 
detector. In Table 1 we have given an estimation of cross sec¬ 
tions for signal and backgrounds considering all possible modes 
with basic cuts. We then investigate the limits on each cou¬ 
pling taking BSM events as the signal. Eor the generation of 
events we use the Monte Carlo event generator MadGraphS [24] 
and the CTEQ6L1 [25] paiton distribution functions. Eurther 
fragmentation and hadronisation are done with a customised 
Pythia-PGS^ [26]. The detector level simulation is performed 
with reasonably chosen parameters using Delphes"* [27] and 
jets were clustered using Fast Jet [28] with the anti-kj- al¬ 
gorithm [29] using the distance parameter, R - 0.4. The fac¬ 
torisation and renormalisation scales for the signal simulation 
are fixed to the Higgs boson mass, m/, = 125 GeV. The back¬ 
ground simulations are done with the default MadGraphS dy¬ 
namic scales. The polarisation is assumed to be -80%. 


*In pp collider like the LHC, the main challenge of this search is to distin¬ 
guish the signal of four bottom quarks in the final state (that hadronise into jets 
(fo-jets)) from the QCD multijet backgrounds. Such challenges and difficulties 
are discussed and performed in ATLAS and CMS studies [20-22]. 

^ We cross checked the modelling of photo-production cross sections from 
MadGraphS by switching on the “Improved Weizsacker-Williams approxima¬ 
tion formula” described in Ref. [23] to give the probability of photons from the 
incoming electron, versus the expectation of the Pythia Monte Carlo genera¬ 
tor. 

^ In Pythia-PGS we modified several parameters in a way to use it for e^p 
collision and to get all required numbers of events demanded in each simula¬ 
tion. The coordinate system is set as for the HERA experiments, i.e. the proton 
direction defines the forward direction. The modifications have been success¬ 
fully validated using neutral current DIS events and switched off QCD ISR. For 
p collisions multiple interactions and pile-up are expected to be very small 
and ai'e switched off in our studies. 

^ For Delphes we used the ATLAS set-up with the modifications in the 
\r]\ ranges for forward and fo-tagged jets up to 7 and 5 respectively with 70% 
tagging efficiency of t)-jets as mentioned in the text. The resolution parameters 
for energy deposits in the calorimeters are based on the ATLAS Run-1 perfor¬ 
mance. 


Process 

cc (fb) 

NC (fb) 

PHOTO (fb) 

Signal: 

2.40 X 10“* 

3.95 X 10-2 

3.30 X 10-“ 

bbbb j: 

8.20 X 10“' 

3.60 X 10+3 

2.85 X 10+3 

bbjjj: 

6.50 X 10+3 

2.50 X 10+“* 

1.94 X 10+“ 

ZZjiZ^bb): 

7.40 X 10“* 

1.65 X 10-2 

1.73 X 10-2 

ttj (hadronic): 

3.30 X 10“* 

1.40 X 10+2 

3.27 X 10+2 

ttj (semi-leptonic): 

1.22 X 10“* 

4.90 X 10+* 

1.05 X 10+2 

hbb j (h —> bh): 

5.20 X 10“* 

1.40 X 10“ 

2.20 X 10-2 

hZj (Z, h —> bb): 

6.80 X 10“* 

9.83 X 10-3 

6.70 X 10-3 


Table 1: Cross sections of signal and backgrounds in charged current 
(cc), neutral current (nc) and photo-production (photo) modes for = 
60 GeV and Ep = 50 TeV, where j is light quarks and gluons. For this 
estimation we use basic cuts \r]\ < 10 for light-jets, leptons and b- 
tagged jets, pr > 10 GeV, AR,ni„ = 0.4 for all particles. And electron 
polarisation is taken to be -0.8. 

3.1. Cut-based optimisation 

We base our simulation on the following kinematic selections 
in order to optimise the significance of the SM signal over all 
the backgrounds; (1) At least four k-tagged jets and one addi¬ 
tional light jet are selected in an event with transverse momenta, 
Pt, greater than 20 GeV. (2) Eor non-fe-tagged jets, the absolute 
value of the rapidity, I 77 I, is taken to be less than 7, whereas 
for Zt-tagged jets it is less than 5. (3) The four k-tagged jets 
must be well separated and the distance between any two jets, 
defined as AR - + (Aifp-, <p being the azimuthal an¬ 

gle, is taken to be greater than 0.7. (4) Charged leptons with 
Pt > \0 GeV are rejected. (5) Eor the largest pj forward 
jet J (the no«-Z 7 -tagged jet after selecting at least four k-jets) 
rfj > 4.0 is required. ( 6 ) The missing transverse energy, $t, 
is taken to be greater than 40 GeV. (7) The azimuthal angle 
between Bj and the Z 7 -tagged jets are: leading jet > 0.4 

and A(l)^^_ sub-leading jet > 0.4. (8) The four / 7 -tagged jets are 
grouped into two pairs such that the distances of each pair to 
the true Higgs mass are minimised. The leading mass contains 
the leading pr-ordered b-']&t. The first pair is required to be 
within 95-125 GeV and the second pair within 90-125 GeV^. 
(9) The invariant mass of all four ^-tagged jets has to be greater 
than 290 GeV. 

In the selections (described above) the Z7-tagging efficiency 
is assumed to be 70%, with fake rates from c-initiated jets and 
light jets to the k-jets of 10% and 1% respectively. Correspond- 


^ Among the four fi-tagged jets, choices of pairing have been performed 
via appropriate selection of mass window, keeping in mind to reconstruct the 
Higgs boson mass, m/j, in the signal as well as the Z-boson mass, mz, in the 
backgrounds. We choose the pair in which the quadratic sum (mi - nic) and 
(m 2 - rric) is smallest, and in each mass m/, mass mi has the largest pj /?-jet, 
rric = (mij - mo) GeV, and normally mo « 20-40 GeV (which is not important, 
since the false pairing will have a much higher quadratic sum). 
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Cuts / Samples 

Signal 

4fo-i-jsls 

2in-jets 

Top 

ZZ 

bbH 

ZH 

Total Bkg 

Significance 

Initial 

2.00 X 10^ 

3.21 X 10’ 

2.32 X lO'’ 

7.42 X 10'’ 

7.70 X 10’ 

1.94 X 10“ 

6.97 X 10’ 

2.36 X 10“ 

0.04 

At least Ab + \ j 

3.11 X 10^ 

7.08 X lO'^ 

2.56 X 10“ 

9.87 X 10’ 

7.00 X 10’ 

6.32 X 10’ 

7.23 X 10’ 

1.08 X 10’ 

0.94 

Lepton rejection pj > 10 GeV 

3.11 X 10^ 

5.95 X lO'^ 

9.94 X 10’ 

6.44 X 10’ 

6.92 X 10’ 

2.26 X 10’ 

7.16 X 10’ 

7.75 X 10“ 

1.12 

Forward jet > 4.0 

233 

13007.30 

2151.15 

307.67 

381.04 

46.82 

503.22 

16397.19 

1.82 

$T > 40 GeV 

155 

963.20 

129.38 

85.81 

342.18 

19.11 

388.25 

1927.93 

3.48 

> 0-4 

133 

439.79 

61.80 

63.99 

287.10 

14.53 

337.14 

1204.35 

3.76 

m‘, 6 [95,125], ml^ e [90,125] 

54.5 

28.69 

5.89 

6.68 

5.14 

1.42 

17.41 

65.23 

6.04 

ma,h > 290 GeV 

49.2 

10.98 

1.74 

2.90 

1.39 

1.21 

11.01 

29.23 

7.51 


Table 2: A summary table of event selections to optimise the signal with respect to the backgrounds in terms of the weights at 10 ab“*. In the 
first column the selection criteria are given as described in the text. The second column contains the weights of the signal process —» hhjve, 

where both the Higgs bosons decay to bb pair. In the next columns the sum of weights of all individual prominent backgrounds in charged current, 
neutral current and photo-production are given with each selection, whereas in the penultimate column all backgrounds’ weights are added. The 
significance is calculated at each stage of the optimised selection criteria using the formula .S = ^2[(S -l-6)log(l + 5/fi) - S], where 5 and B 
are the expected signal and background yields at a luminosity of 10 ab“*respectively. This optimisation has been performed for Eg = 60 GeV and 
E„ = 50 TeV. 




Figure 2: Azimuthal angle distributions, at Delphes detector-level, between missing transverse energy, tj, and the forward jet, J, in the SM 
(including backgrounds) and with the anomalous hhh, hWW and hhWW couplings. The error bars are statistical. 


ing weights® at a particular luminosity of 10 ab^'for a signal, 
and all backgrounds with significance has been tabulated in Ta¬ 
ble 2. Significance at all stages of the cuts are calculated using 
the Poisson formula S - -\/2[(S + B) log(l + SIB) - 5], where 
S and B are the expected signal and background yields at a par¬ 
ticular luminosity respectively. 

From Table 2 we recognise that selection on the forward jet 
in the FCC-he type machine plays a very significant role in 
distinguishing the signal with respect to background. By se¬ 
lecting events with tjj > 4.0, there is loss of 25% of signal 
events, while the total background loss is around 80%. The 
next significant cut on missing energy {$t >40 GeV) is also 
very significant as due to this cut there is loss of 88% of events 


®Here weights mean the number of expected events at a particular luminos¬ 
ity. The number of events of the photo-production of 4h-t-jets are derived using 
the efficiencies of the Monte Carlo samples due to the low statistics. The other 
backgrounds are obtained directly from the event selection. 


in total background, however, for the signal there is a loss of 
only 34% of events after forward jet selection. Furthermore the 
mass window cut for the invariant mass of two fi-tagged jets, 
after the > 0.4 selection, significantly reduces the total 

background events to 5%, only while the signal events remains 
at 40%. Efficient requirements on the invariant mass window 
of four fi-jets are efficient, such that to reduce backgrounds by 
44% leads to a signal of 90% with respect to the previous two 
fi-tagged jet mass window selection. And hence there is a 20% 
enhancement in the significance obtained in comparison to the 
two different mass window selection criteria, and overall with 
respect to initial events this cut-based optimisation is enhanced 
from a 0.04 to 7.51 significance. Here it is also important to 
mention that photo-production of the Ab final state is one of the 
main background with similar topological final states from Zh, 
where Z, —> bb, and is equally important. Hence choice of 
efficient selection criteria is too important to reduce these back¬ 
grounds. 
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3.2. Kinematic distributions and observable 


For our analysis we take ad hoc values of positive and neg¬ 
ative couplings in such a manner that the production cross sec¬ 
tion does not deviate much from the SM value, and in particu¬ 
lar modifications in the shapes of the normalised azimuthal an¬ 
gle distribution between the missing transverse energy and the 
leading (forward) jet are studied, in addition to other kinematic 
distributions. 

Taking into account all the above criteria we study 
BSM modifications in various differential distributions at the 
Delphes detector-level. This leads to the following observa¬ 
tions: ( 1 ) pt has the usual tail behaviour, i.e. the number of 
events are more populated in the higher pj region with respect 
to the SM for the chosen values of the anomalous couplings. 
(2) In cases of the rf distributions: (a) For the forward jet, par¬ 
ticularly for the couplings of hWW and hhWW vertices, the 
mean rj is more central in the detector. The behaviour is similar 
if we increase the c.m.s. energy of the collider by increasing 
to higher (> 60 GeV) values. For hhh couplings the rj distri¬ 
bution remains the same as for the SM. (b) In case of /j-tagged 
jets, for all values of anomalous couplings, the distribution is 
populated around the value of 77 of the SM distribution. (3) For 
the specific observable of the azimuthal angle difference be¬ 
tween missing transverse energy and the forward jet (A(pfj.j) 
the shapes are clearly distinguishable from the SM. This be¬ 
haviour is shown in Fig. 2, where the values of the couplings 
are ad hoc. However, these values are taken only for the pur¬ 
pose of illustration, and in the limit of the couplings going to 
their SM values the shapes will coincide with the SM distri¬ 
butions. The specific characteristics of the curves also depend 
on the details of the selection requirements, but the qualitative 
differences could be seen at every selection step. The shape of 
the curves is due to the fact that all new physics couplings have 
a momentum dependent structure (apart from and positive 
or negative interference with SM events. We note that A(pgj.j 
is a novel observable and commands more focused and deeper 
analyses. In this regard one should follow the analysis (as per¬ 
formed in Ref [30]) based on an asymmetry with two equal 
bins in ^ njl, defined as 

_ |AA0>7r/2l - \A^lp<n:|2\ 

|Aa^>./2| + |Aa^<./2|’ 

in which A is the yields obtained for the given ah 'data af¬ 
ter all of the selections, including both the signal and back¬ 
grounds. Table 3 shows the estimation of the asymmetry for 
a set of representative values of the couplings, shown in Fig. 2, 
along with the associated statistical uncertainty. Though the 
new physics couplings are representative, we can infer the sen¬ 
sitivities of these couplings with respect to the SM-i-Bkg esti¬ 
mation of asymmetry from Table 3, where seems to have 
large fluctuations for both positive and negative choices of its 
values. Similarly the sensitivities of g,® and can also be 
noted, however g'*' is more sensible for the negative choice of 

‘-'/ifiVvvv ^ 

its value. The study of the sensitivity of the non-standard cou¬ 
plings through this asymmetry observable considering a kine¬ 
matic distribution is basically corresponding to two bins by 


Samples 


o-(fb) 

SM-bBkg 

0.277 + 0.088 


„(i) 

^hhh 

= 1.5 

= 2.0 

0.279 + 0.052 

0.350 + 0.053 

0.18 

0.21 

„( 2 ) 

°hhh 

II II 

p 0 

0.381 +0.050 

0.274 + 0.024 

0.19 

0.74 

„(i) 

^hWW 

II II 

p 0 

0.506 + 0.022 

0.493 + 0.020 

0.88 

0.94 

„( 2 ) 

= - 0.02 

= 0.02 

0.257 + 0.025 

0.399 + 0.040 

0.67 

0.33 

8hWW 

p 0 

I 

II II 

0.219 + 0.016 

0.228 + 0.016 

1.53 

1.53 

„(i) 

°hhWW 

= - 0.05 

= 0.05 

0.450 + 0.033 

0.254 + 0.029 

0.52 

0.68 

„( 2 ) 

^hhWW 

= - 0.03 

= 0.03 

0.462 + 0.022 

0.333 +0.018 

1.22 

1.46 

ShhWW 

II II 

p p 

0.351 +0.020 

0.345 + 0.020 

1.60 

1.61 


Table 3: Estimation of the asymmetry, defined in Eq. (9), and statis¬ 
tical error associated with the kinematic distributions in Eig. 2 at an 
integrated luminosity of 10 ab“'. The cross section (cr) for the corre¬ 
sponding coupling choice is given in the last column with same param¬ 
eters as in Table 1 . 

dividing the whole distribution in two halves with large bin- 
width. Moreover, this kind of study can be further appended 
with finer bin widths and a more efficientanalysis (for exam¬ 
ple in Ref. [30]). However, these detailed analyses are beyond 
the scope of this article. 

3.3. Exclusion limits through fiducial cross section as a func¬ 
tion of luminosity 

Furthermore, we probe the exclusion limits on these cou¬ 
plings as a function of the integrated luminosity, with the log- 
likelihood method described in Ref. [31], using directly the 
fiducial inclusive cross section as an observable. In Fig. 3 we 
present exclusion plots at 95% C.L. for anomalous hhh, hWW 
and hhWW couplings, where the shaded areas are the allowed 
region. The exclusion limits are based on the SM ‘di-Higgs sig¬ 
nal -H backgrounds’ hypotheses considering BSM contributions 
as the signal at the given luminosity. Each limit is given by 
scanning one coupling and fixing the other couplings to their 
SM value, where a 5% systematic uncertainty is taken into ac¬ 
count on the signal and background yields respectively. From 
Fig. 3 our observations are as follows: (1) If the integrated lu¬ 
minosity exceeds 0.5 ab“' gjjj is restricted to be positive, 
is allowed to be within 0.7-2.5 when the integrated luminosity 
reaches 15 ab“* as for values of 1 < g'*' <2.1 the cross sec- 
tion is smaller than the SM di-Higgs production. (2) The g® 
value is restricted to around 10 *. We only exclude the posi¬ 
tive part of this coupling because its negative part has cancel¬ 
lations with the SM di-Higgs cross section. (3) The sensitivity 
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Figure 3: The exclusion limits on the anomalous hhh (top panel), hWW 
(middle panel) and hhWW (lower panel) couplings at 95% C.L. as a 
function of integrated luminosity (shaded areas). Note that the allowed 
values of g® and are multiplied by 5 and 10 respectively to high¬ 
light their exclusion region, since the values are of the order 10“'. 

for hWW couplings, namely and can be better probed 
at much lower energies and luminosity at the LHeC using the 
single Higgs boson production as shown in Ref. [18]. However, 
we have shown the sensitivity of which is not considered 
in Ref. [18], and it is of the order 10“^ in the allowed region. 
(4) One important aspect of di-Higgs production in this type of 
collider is that one can measure the sensitivity of the hhWW 
couplings also. In our analysis, since the CP-even (odd) cou¬ 
pling (8i,i,ww') has similar Lorentz structures, with the sen¬ 
sitivity of the exclusion plot having almost the same order of 
magnitude. However, the structure of allows a compar¬ 
atively narrower region of values. The couplings belonging to 
both the hWW and hhWW vertices are strongly constrained be¬ 
cause of their high production cross section at very low values 
of the couplings. By increasing the luminosity from 0.1-1 ab“' 
the constraint on the couplings increases and its limits are re¬ 
duced by a factor two. A further increase of the luminosity 
will not change the results. All limits are derived by varying 
only one coupling at a time, as mentioned earlier. The exclu¬ 
sion limits on the couplings in this analysis are based on the 
constraints from an excess above the SM expectation while po¬ 
tential deficits from interference contributions are not sensitive 
yet to be used for limit settings. 

3.4. Prospects at higher Eg and sensitivity of the Higgs self 
coupling 

Finally we discuss what happens once the electron energy 
is increased to higher values, where we focus our analysis on 


a determination of the SM Higgs self coupling, assuming no 
further BSM contributions. Without going into detail we can 
note that with increasing E^ (from 60 GeV to 120 GeV) the SM 
signal and dominant background production cross sections are 
enhanced by a factor of 2.2 and 2.1 respectively. As a result, the 
cut efficiency for the selection of four Z7-tagged jets and one for¬ 
ward jet is improved, but for the other cuts described previously 
(invariant mass, $t, W and ^<p$jf) it remains very similar. This 
leads to an enhancement of the selected signal and dominant 
background events by a factor 2.5 and 2.6 respectively. Hence 
we would obtain the same statistical precision with only 40% 
of the luminosity of an Eg - 60 GeV beam when increasing 
the electron energy to 120 GeV. At an ultimate integrated lu¬ 
minosity of 10 ab“', increasing Eg from 60 to 120 GeV would 
increase the significance of the observed SM di-Higgs events 
from 7.2 to 10.6, obtained from a log-likelihood fit. This in¬ 
cludes a 5% signal and background systematics mentioned ear¬ 
lier. For the SM Higgs boson self coupling, where the scal¬ 
ing factor is expected to be gjjj^ = 1, we perform an intelligent 
signal injection test, which gives locally measured uncertain¬ 
ties for From this test the Icr error band around the ex- 

pected SM strength of this coupling is ^-O^-o i 7 (o i 2 ) 

Eg = 60(120) GeV. 

4. Validity of EFT 

In the EFT-based approach for our analyses, the usual SM 
Lagrangian is supplemented by higher-dimensional operators 
that parametrise the possible effects of new states assumed to 
appear at energies larger than the effective scales identified with 
niw (or equivalently v) by restricting the operators of dimen¬ 
sion less than or equal to six. We have estimated the sensitivity 
of the involved coupling coefficients appearing in the effective 
Lagrangians in Eqs. (2)-(4) with the EW scale for the deriva¬ 
tive terms. A detailed discussion with general couplings and 
mass scales with a higher dimension EET Lagrangian can be 
found in Ref. [32]. Eor the processes at high energy, it is well 
known that an EET approach provides an accurate description 
of the underlying new physics as long as the energies are be¬ 
low the new physics scale. A, and thus the limits on the cou¬ 
plings obtained in the above analyses shall degrade at scales 
higher than the EW scales (since for the fixed values of the 
couplings, the interference and pure BSM terms always give 
low contributions in the cross section measurements for high 
values of scale choice). Also for 0(1) values of anomalous cou¬ 
plings (apart from §(■■■) TeV-scale momenta, one 

reaches the regime where the operators in Eqs. (2)-(4) may not 
be dominant, and operators with four and more derivatives may 
be equally important. In other words, the EFT behind these 
Lagrangian expansions breaks down. It would be important 
then to know how much the projected sensitivity depends on 
events that violate this EFT bound. With an EW precision test 
in Ref. [33], it is shown how an EET’s reach deteriorates when 
only data below the cutoff scales are employed on the mass vari¬ 
able in the case of Drell-Yan processes at the LHC. 

A similar exercise can be performed in our case to estimate 
the deterioration of limits on anomalous tensorial couplings 
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Figure 4: Percentage of deterioration of exclusion limits of anomalous 
tensorial couplings (shown in Fig. 3) with respect to the upper di- 
Higgs invariant mass cut = rrf^ [in GeV] for fixed luminosity of 
1 ab^'{blue) and 10 ab“‘(rerf). The numbers in the vertical axis above 
(below) 0 is the degradation in the upper (lower) limits. 

.g® j and g(...) (the coupling coefficients which corresponds to 
momentum-dependent Lorentz structure) as a function of the 
cut-off scale. In this approach we put an upper cut on the 
di-Higgs invariant mass such that EFT-violating events 

(> m^^) are cut away, and then we estimate by how much 
the projected sensitivity of g®, and g„,,„,, g, 

degrades with respect to their previous limits. In Fig. 4 we 
present the percentage of deterioration of the exclusion lim¬ 
its of these anomalous effective couplings by selecting events 
below e [0.35,1] TeV for fixed luminosity of 1 ab^'and 
10 ab"^at 95% C.F. It is apparent from Fig. 4 that the deteri¬ 
oration in the limits of these anomalous couplings is large for 
low values of the cut, because the effective cross section 
decreases (which is equivalent to the increase of the scale A 
of the tensorial couplings) with the decrease of the values of 
Comparing the exclusion limits obtained in Fig. 3 we ob¬ 
served that at m'^i^ - 350 GeV the percentage of deterioration 
in g*^* is more than 100%, while other hWW and hhWW cou- 
plings deteriorate by 60-80% on both upper and lower sides at 
1 ab“*and 10 ah '. After 350 GeV a sudden decrease in degra¬ 
dation percentage can be noticed for = 400 - 450 GeV for 
all couplings. Furthermore around 500 GeV for g®, it remains 


^Note that in previous subsections we used the notation nm, for the lowest 
cut on di-Higgs invariant mass. Here we use to avoid confusion since for all 
analyses apart from the EFT validity we selected the events maj > 290 GeV to 
suppress backgrounds and increase the overall significance. Here, to investigate 
the sensitivity of BSM tensorial couplings, we chose the events below cuts, 
keeping maj > 290 GeV so that an one to one comparison can be performed. 


18% while others are around 10%. Beyond a 650 GeV cut, all 
the couplings converge to the original value of limits obtained 
in our previous analyses, as shown in Fig. 3. 

5. Summary and conclusions 

We conclude that the FCC-he, with an FRF energy of Ee > 
60 GeV and a proton energy Ep - 50 TeV, would provide sig¬ 
nificant di-Higgs event rates, and through this channel one can 
probe accurately the Higgs boson self coupling provided that 
integrated luminosities of more than 1 ab“^ may be achieved. 
Along with the Higgs self coupling one can search for any 
BSM signal through the measurement of the anomalous hhWW 
contributions. One interesting feature of this type of machine 
is recognised by identifying forward jets in the signal events 
where an appropriate selection, as shown for our study, reduces 
backgrounds efficiently around 80% with a loss of only 25% of 
signal events. Our work also shows that A(pgj.j is a very good 
observable for any new physics contributions in the given chan¬ 
nel. Estimation of an asymmetry observable in for this 

kinematic distribution gives a preliminary idea of sensitivities 
of any new non-standard couplings. The limits on each cou¬ 
pling are set by measuring the observed event rate. But the 
asymmetry in can provide more distinguishability of the 

new physics, especially cancelling many potential systematics, 
which is helpful to distinguish the signatures of each model. An 
exclusion limit with respect to luminosity for these couplings is 
studied, and a signal injection test shows the uncertainty of the 
Higgs self coupling around its expected SM strength. 

With all these analyses we infer that the order of sensitiv¬ 
ities of all non-standard couplings considered for our study 
within most of the luminosity ranges are consistent with the 
adopted methodology of asymmetry observable, and exclusion 
limits through fiducial cross sections at 95% C.F. However, at 
luminosities ~10-15 ah 'or higher, the method based on fidu¬ 
cial cross sections constrains the non-standard couplings more 
tightly. In addition to the fiducial cross sections, the drastic 
change in shape for around 2.0 with respect to the 

SM in Fig. 2 (left), suggest that using further observables like 
may significantly improve the sensitivity of BSM cou¬ 
plings in general. It is to be noted that the non-standard momen¬ 
tum dependent structures of the EFT breaks down at the TeV 
energy regime for couplings of (9(1) and then additional deriva¬ 
tive terms become relevant. Hence we also show the deteriora¬ 
tions in the limits of anomalous tensorial couplings for different 
regions of di-Higgs invariant mass (upper) cuts by an exclusion 
at fixed luminosity of 1 ah 'and 10 ab“'corresponding to 95% 
C.F. limits, with respect to the limits obtained using fiducial in¬ 
clusive cross section as an observable. This method is used as 
an alternative approach to estimate the sensitivity of the scale 
dependent couplings in EFT and gives a probe to understand 
the regions where validity of EFT breaks down. 

Our studies show a unique capability and potential of FCC- 
he collider to probe the precision measurement of not only the 
Higgs boson self coupling but also other involved couplings 
with tensorial structure through di-Higgs boson production. 
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